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TECHNICAL  DISCUSSION 

I.  High  Gain  Lasers 

During  the  year  1972  we  concluded  our  work  on  high  gain  lasers. 

The  work  consisted  of  theoretical  and  experimental  investigation 
of  spectral,  temporal  and  modal  properties  of  extremely  high  gain  lasers. 

The  experiments  were  conducted  using  xenon  3.51  pm  lasers  with  gains  up 
to  60  db/m. 

A  number  of  significant  effects  were  predicted  by  the  theory  and 
checked  in  experiments.  A  number  of  resulting  publications  (see  References 
1  and  2)  the  the  Ph.D.  thesis  of  Dr.  L.  Casperson  summarize  the  results.^ 

II.  Nonlinear  Optics 

Accurate  measurements  were  performed  to  determine  the  nonlinear 
optical  properties  of  a  number  of  important  materials.  These  include  KDP, 
RDP,  RDA,  and  Li  1 03 .  The  results  are  described  in  the  enclosed  Reference  3. 

A  theoretical  analysis  was  carried  out  which  points  out  a  new 
technique  for  phase  matching  in  nonlinear  optical  interactions.  The 
technique  involves  fine  corrugation  of  thin  dielectric  waveguides.  Details 
are  presented  in  Reference  4. 

III.  Optical  Damage 

Working  with  C.  Giuliano  of  the  Hughes  Research  Laboratories,  it 
was  shown  experimentally  and  theoretically  that  the  use  of  an  elliptical 
cross  section  Gaussian  beam  can  lead  to  large  increases  in  the  self  focus¬ 
ing  damage  threshold  of  laser  materials. 
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This  suggests  the  incorporation  of  astigmatic  optics  in  high 


power  laser  systems  (see  Reference  5). 
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Gain  and  Dispersion  Focusing  in  a  High  Gain  Laser 


Lee  W.  Casperson  and  Amnon  Yariv 


The  transverse  modes  of  a  laser  resonator  containing  a  medium  with  a  strong  radial  gain  profile  differ 
greatly  from  the  modes  of  a  similar  resonator  containing  a  low  gain  medium.  Foeusiug  and  defocusmg 
effects  result  from  the  gain  profile  and  from  the  associated  dispersion  profile.  The  dispersion  foeusing 
causes  an  asymmetry  in  the  power  output  as  the  laser  is  tuned  aeross  the  gain  line.  The  theory  has  been 
verified  using  a  high  gain  3.51-p  xenon  laser. 


I.  Introduction 

fr  lenslike  media  having  an  approximately  quadratic 
radial  variation  of  either  the  gain  or  the  index  of  refrac¬ 
tion  there  is  a  focusing  or  defoeusing  of  propagating 
Gaussian  beams.1,2  A  positive  index  profile  results  in  a 
spot  size  of  an  unmatched  beam  which  oseillates  periodi¬ 
cally  with  distance,  while  a  positive  gain  profile  causes 
the  beam  to  approach  a  steady-state  spot  size.  In  a 
high  gain  gas  discharge  there  is  often  a  strong  profile 
of  the  gain,  so  that  gain  focusing  is  likely  to  be  impor¬ 
tant.  However,  whenever  a  medium  has  a  gain  spec¬ 
trum,  it  must  also  have  an  associated  dispersion  spec¬ 
trum.  Therefore,  dispersion  foeusing  effects  must  al¬ 
ways  aeeompany  gain  focusing.  It  is  the  purpose  of 
this  paper  to  show  that  the  two  types  of  focusing  are 
comparable  in  importance  and  that,  moreover,  the  dis¬ 
persion  foeusing  may  lead  to  asymmetry  in  Lamb  dip 
measurements.  Ordinarily  asymmetries  in  the  output 
of  simple  gas  lasers  are  attributed  to  collisions  between 
the  atoms.3,4  The  lasers  considered  are  assumed  to  be 
operated  very  near  threshold,  so  that  self-foeusing  is 
unimportant.6  This  treatment  is  based  on  one  given 
elsewhere.6 

s 

II.  Theory 

A  lenslike  medium  is  one  in  whieh  the  eomplex  propa¬ 
gation  constant  k  varies  quadratically  with  radius  ac¬ 
cording  to 

k  -  kt  —  ]ktr *.  (1) 

The  real  and  imaginary  parts  of  k  are  given  by 

k  =  0  +  tor.  (2) 


Both  authors  were  with  the  Division  of  Engineering  mid  Ap¬ 
plied  Science,  California  Institute  of  Technology,  Pasadena, 
California  91109;  L.  W.  Casperson  is  now  with  the  UCLA 
School  of  Engineering  and  Applied  Science. 
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Here  a  is  the  exponential  gain  constant  for  the  electrie 
field  and  /S  is  related  to  the  index  of  refraetion  by 

0  =  2*n/X.  (3) 

The  easiest  way  to  study  the  propagation  of  Gaussian 
beams  through  such  lenslike  media  is  in  terms  of  the 
complex  beam  parameter  q  given  by 

1/9  =  (l/tf)  -  t(\m /**>>*),  G) 

where  R  is  the  radius  of  eurvature  of  the  phase  fronts, 
w  is  the  spot  size  of  the  beam,  and  Am  is  the  wavelength 
in  the  medium.  The  same  beam  parameter  may  also 
be  used  to  characterize  higher  order  Laguerre-Gaussian 
and  Hermite-Gaussian  beams. 

From  the  wave  equation  one  finds  that  the  propaga¬ 
tion  of  q  is  governed  by  the  Ricatti  equation1 

(1/9)*  +  (d/dz)(l/9)  +  (V*.)  =  0-  (5) 

For  simplicity  we  assume  here  that  saturation  is  unim¬ 
portant  and  the  coefficients  ko  and  ki  are  independent 
of  position  along  the  axis  of  the  medium.  The  solution 
of  Eq.  (5)  is 

l/9i=  U(l/9i)  cos(Jt*/Jt«)**l  -  ((*»/*«)*  sin(fci/fc0)**] } 

4-  { 1(1  /9i)(fco/fr»)*  sin(fc»/fco)**l  +  (cos(fci/fco)L]}.  (6) 

Exeept  when  (k«/k0)'  is  real,  the  beam  parameter  at 
large  distances  approaches  the  limiting  value 

+  {l/9»  =  TiffciA'o)  Mi 

where  the  upper  sign  is  used  if  Im(fcj/A>)*  >  0  and  the 
lower  sign  is  used  if  Im(t:/to)*  <  0- 
In  the  special  ease  Im(/.';/A'o)  *  =  0,  Re(Aj/Ao)  7^  0 
(no  gain  profile),  the  beam  parameter  oseillates  periodi¬ 
cally  without  approaching  a  steady-state  value.  If 
(kt/kt)1  =  0  so  that  there  is  no  profile  at  all,  Eq.  (6) 
reduees  to  the  free  space  result 

91  =  91  +  z-  W 

Otherwise,  after  some  oseillation  the  beam  goes  to  the 
limit  given  in  Eq.  (7).  From  Eqs.  (4)  and  (7)  the  limit- 


462  APPLIED  OPTICS  /  Vol.  11,  No.  2  /  February  1972 


mg  spot  size  is 

»•  =  (X„/5r)*[±Re(A*/Ar0)*]  —  *  (9) 

and  the  limiting  radius  of  curvature  of  the  phase  fronts 
is 

Rm  **  (iliu^/M1!"'.  (10) 

The  spot  size  given  by  Eq.  (9)  must  be  real  if  the 
beam  is  to  be  confined  by  the  medium.  With  our  sign 
convention  this  means  that  the  real  and  imaginary  parts 
of  (A-2/A-0) 1  must  have  the  same  sign,  or  that 

Im (ki/k„)  >  0.  (11) 

With  Eq.  (2)  this  stability  condition  is  simply 

«« >  0,  (12) 

provided  that  the  gain  per  wavelength  is  small  (ft,  » 
oro)-  In  other  words,  the  gain  profile  causes  a  damped 
focusing  such  that  the  beam  approaches  a  stable  steady- 
state!  value  if  and  only  if  the  gain  is  highest  on  the  axis 
of  the  medium.  It  is  only  this  stable  situation  which  is 
of  interest  in  the  remainder  of  this  paper.  From  Eqs. 
(9)  and  (10)  the  spot  size  is 

=  (A„/ir)J{  ±Re[(0i  -f  lai)/0ol*l  -1  (13) 

and  the  radius  of  curvature  of  the  phase  fronts  is 

ft.  =  {±lm[(0,  +  «'ari)//SoI  1) (14) 

We  have  described  so  far  the  propagation  of  Gaussian 
beams  in  long  high  gain  laser  amplifiers.  The  transverse 
inodes  of  a  laser  resonator  containing  a  high  gain  me¬ 
dium  are  most  easily  found  using  complex  beam 
matrices.2  In  the  remainder  of  this  paper  we  consider 
only  the  simplest  type  of  high  gain  laser,  which  consists 
of  a  plane  parallel  resonator  filled  with  a  lenslike 
medium.  Such  a  configuration  can  be  well  approxi¬ 
mated  in  practice.  The  spot  size  everywhere  in  this 
laser  is  equal  to  the  limiting  value  given  by  Eq.  (13), 
which  can  be  written  as 

w  ~  ((*,“«/4A.){  (1  +  (0, /a,)*]*  +  (0j/<*j)j)_1.  (15) 

This  result  can  be  applied  to  lasers  with  various  types 
of  index  profiles.  In  the  following,  only  profiles  due  to 
the  gain  and  dispersion  will  be  considered. 

To  find  how  the  spot  size  in  Eq.  (15)  depends  on  fre¬ 
quency,  it  is  necessary  to  know  the  frequency  depen¬ 
dence  of  a2  and  ft>.  In  a  Doppler-broadened  medium 
the  unsaturated  intensity  gain  g  is  a  Gaussian  function 
of  frequency  and  a2(x)  may  be  written 

«>(*)  =  (0172)6-*’,  (16) 

where  x  =  2(v  —  i/0)(ln2),/Ai'D  is  a  normalized  fre¬ 
quency  and  g-i  is  the  line  center  value  of  the  quadratic 
term  in  the  intensity  gain  constant.  The  quadratic 
term  in  the  index  of  refraction  of  a  Doppler-broadened 
medium  is7 

«i(x)  =  (\g,'/2*^)F(x),  (17) 

where  F(x)  is  Dawson’s  integral  given  by 

F(x)  =  «-*’  J  e‘'dt.  (18) 


Therefore  ft, (a;)  is  simply 

/»•(*)  =  (9i'/^)F(x).  (19) 

Using  Eqs.  (16)  and  (19)  in  Eq.  (15),  one  obtains  the 
following  expression  for  the  frequency-dependent  spot 
size  of  a  plane  parallel  gain-focused  laser: 

—[£'*•■■  (i- mr+?eF')]" 

(20) 

It  is  convenient  to  define  the  normalized  spot  size  iu*(x) 
given  by 


. . ((,+[^]')7^n  *> 


w'hich  is  equal  to  the  actual  spot  size  divided  by  its 
line  center  value.  A  plot  of  w*(x)  is  given  in  Fig.  1. 


Fig.  1 


Also  plotted  in  the  figure  is  the  function  eI,/4,  which 
would  represent  the  frequency  dependence  of  the 
normalized  spot  size  if  dispersion  focusing  were  ne¬ 
glected.  The  asymmetry  of  the  spot  size  results  from 
the  positive  and  negative  dispersion  focusing,  which 
occur,  respectively,  for  positive  and  negative  values  of 
the  frequency  x.  The  minimum  value  of  the  spot  size 
is  at  a  frequency  of  about  x  ~  0.6  rather  than  at  line 
center. 

We  are  primarily  interested  here  in  inhomogeneously 
broadened  media,  but  the  results  for  homogeneous 
broadening  are  similar.  In  an  unsaturated  homoge¬ 
neous  medium  the  gain  is  the  Lorentzian 


«» (v)  -  (0»72){i/n  +  »■)!,  (22) 

wheie  y  —  2(v  —  va)/bvn  is  a  normalized  frequency. 
The  index  of  refraction  is 


My)  =  (Aff»74»  )ls//(l  +  !/*)].  (23) 

Therefore,  from  Eq.  (15)  the  spot  size  is  given  by 

u>  =  |(it0,78a»)(1/(1  +  !/*)Jl(l  +  *•)*  +  V]J  (24) 
These  homogeneous  results  are  not  considered  further. 


-15  -1.0  -.5  0  .5  1.0  1.5  x 

Fig.  1.  Solid  line  is  the  normalized  spot  size  as  a  function  of 
frequency.  Dashed  line  is  the  spot  size  neglecting  dispersion 
focusing. 
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The  difference  between  the  spot  sizes  indicated  by 
the  two  curves  in  Fig.  1  is  at  most  of  the  order  of  10%. 
Nevertheless,  the  dispersion  focusing  eould,  in  principle, 
bi  deli  £lpd  directly  by  scanning  the  output  beam  |_*rnfile 
of  a  laser  if  the  oscillation  frequency  were  known.  A 
more  important  practical  consequence  of  the  dispersion 
lueuui^i  i  mi  rtfc  a  rant  piuarr  juaiaw  .f  bam 
becomes  asymmetric  about  line  center.  One  expects 
that  for  equal  values  of  the  gain,  the  larger  the  beam 
spot  size,  the  greater  the  output  power.  In  particular, 
the  maximum  power  output  of  the  lasc  does  not  occur 
when  the  oscillation  frequency  is  at  gain  center  (neglect¬ 
ing  the  Lamb  dip,  of  course).  The  asymmetry  in  the 
power  output  is  easily  measured  and  provides  a  fairly 
direct  indication  of  dispersion  foeusing. 

The  growth  of  intensity  in  a  saturating  one-dimen¬ 
sional  inhomogeneously  broadened  laser  is  governed  by 
the  well-known  expression 

dl/dz  =  gI/(  1  +  sl)>,  (25) 

where  s  is  a  saturation  parameter.  For  simplicity  dis¬ 
tributed  losses  are  neglected.  For  very  weak  satura¬ 
tion  (s I  «  1)  Eq.  (25)  reduces  to  the  homogeneous  ap¬ 
proximation 

dl/dz  =  gI/{  1  +  Is/).  (26) 

If  the  laser  beam  were  uniform  over  an  area  of  radius  w 
and  there  were  no  gain  profile,  then  the  corresponding 
result  for  the  growth  of  the  total  power  would  be 

dP/dz  =  g„P/ll  +  KsP/irui1)!.  (271 

Equation  (27)  can  be  shown  to  be  also  valid  for  a 
Gaussian  beam  with  nearly  plane  phase  fronts  in  a 
medium  with  a  quadratic  gain  profile,  provided  that 
the  square  of  the  spot  size  iv  is  much  smaller  than  the 
square  of  the  discharge  diameter  r0.6  This  is  not  always 
a  good  approximation  but  it  is  valid  for  our  experi¬ 
ments.  0o  is  the  unsaturated  gain  at  the  axis  of  the 
medium.  The  result  for  a  homogeneously  broadened 
medium  would  be  the  same  as  Eq.  (27)  except  that  the 
factor  £  would  be  missing. 

If  one  mirror  is  highly  reflecting,  Eq.  (27)  can  be  in¬ 
tegrated  for  one  loop  through  the  laser  medium,  and  the 
result  is 

1(»/Tt «’)(/',  -  P.)  =  2gd  -  In (P,/P,),  (28) 

where  l  is  the  length  of  the  medium.  This  integration  is 
possible  as  long  as  the  right  and  left  traveling  beams 
interact  ',\i"t1l  different  velucn.y  Classes  uf  illOiilfc.  If  the 
reflectivity  of  the  output  mirror  is  R  and  the  transmis¬ 
sion  is  T,  then  Eq.  (2S)  may  be  solved  for  the  output 
power  P0  as 

P,  =  (2xt1>V*)[r/(l  -  R)](2gol  +  ln«).  (29) 

A  similar  result  can  De  oDtamea  tor  nomogeneousiy 
broadened  lasers. 

Using  Eq.  (20)  for  the  frequency-dependent  spot 
size  with  0o  =  0t ,'c~x\  one  obtains  finally 

P,(t)  =  (4ro/s)(irXn/1.44(?o')h*,/!(7,/(l  -  ft)l (20,7c-'’  +  In R) 
X  ({1  +  [2P(x)e*,/iril,J »  +  (2P(x)e*Vr*l)_l.  (30) 


We  have  assumed  here  that  the  medium  has  a  Bessel 
function  radial  gain  profile,  so  that  the  quadratic  term 
iss 

01  ”  0o\,2.»»/  To*;.  (31) 

Equation  (30)  is  the  general  expression  for  the  fre- 
of  itn>  anr|m  pwr  dl  i 

broadened  laser  oscillator  (neglecting  the  Lamb  dip). 
The  exact  calculation  of  the  Lamb  dip  is  difficult  in 
high  gain  lasers.  However,  as  long  as  the  homogeneous 
line  width  is  small  compared  to  the  Doppler  line  width, 
the  Lamb  dip  provides  a  useful  indication  of  line  center 
( x  =  0)  without  significantly  affecting  the  over-all  line 
shape. 

It  is  convenient  to  define  the  normalized  power  spec¬ 
trum 

P,*(x)  =  (e-*’  -  6)c**/I({  1  +  [2P(x)e*7*1l,J1 

+  [2P(x)c*7irJ])-‘,  (32) 

where 

6  -  —  (ln/f/20,,7)  (33) 

is  a  threshold  parameter  which  is  less  than  unity  if  the 
laser  is  saturated.  Equation  (32)  is  ploticd  in  Fig.  2 
for  various  values  of  the  parameter  b.  Evidently  when 
the  laser  is  above  threshold,  the  greatest  output  occurs 
at  a  slightly  negative  frequency  rather  than  at  line 
center.  Near  threshold  ( b  -*■  1)  this  effect  diminishes 
and  the  greatest  output  occurs  near  x  =  0.  Also  plotted 
in  Fig.  2  is  the  gain  spectrum  e~z\ 

The  asymmetry  of  the  power  spectrum  is  most  con¬ 
veniently  characterized  by  the  location  in  frequency  of 
the  power  maximum.  In  Fig.  5  is  a  plot  of  this  fre¬ 
quency  as  a  function  of  the  threshold  parameter  b. 
The  plots  in  Figs.  2  and  5  are  not  quantitatively  correct 
for  all  values  of  b.  From  Eqs.  (27)  and  (29)  this  homo¬ 
geneous  approximation  is  only  valid  in  an  inhomoge- 
ncously  broadened  medium  as  long  as  the  product  (2 g0l 
+  lni?)(l  —  R)~l  is  small  compared  to  unity.  In  a 


Fig.  2.  Normalized  power  output  as  a  function  of  frequency  for 
various  values  of  the  threshold  parameter.  Dashed  line  is  the  gam 
spectrum. 
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high  loss  laser  (R  <J£  1)  this  is  the  same  as  requiring  that 
b  be  nearly  equal  to  unity.  Nevertheless,  these  results 
are  expected  to  be  qualitatively  correct  for  most  values 
of  6. 

It  is  possible  to  obtain  a  simpler  approximate  expres¬ 
sion  for  the  output  power  as  a  function  of  frequency 
which  is  valid  for  small  values  of  x.  One  finds  after 
some  algebra  that  Eq.  (32)  may  be  expanded  to  second 
order  in  a;  as 


P,*(x)  ~  (i  -  6)  _  [(i  _  b)/ri]x  +  { [(1  _  6)/2t] 

-  1(1  +  6)/2J}i*.  (34) 
iht  iirtHwnWki  oi  this  sptetruin  occurs  ax  the  irequency 
=  {r~1  ~  1(1  +  b)/(l  -  6)]xiJ-i.  (35) 

If  i  is  nearl  k  a^jal  to  unity  llu  o  0*  fe  ^ 


Xmmx - [(1  —  6)/2t*].  (36) 

This  equation  is  plotted  as  a  dashed  line  in  Fig.  5.  The 
two  line.'  in  the  figure  are  in  good  agreement  hi  the  limit 
of  weak  saturation. 

In  this  section  we  have  discussed  focusing  effects 
which  are  due  to  the  gain  and  dispersion  profiles  which 
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Fig.  3.  Experimental  setup. 
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may  be  associated  with  a  high  gain  laser  transition.  It 
vas  shown  that  the  spot  size  of  a  waveguided  beam  is 
greater  for  frequencies  below  gain  center  than  for  fre¬ 
quencies  greater  than  gain  center.  As  a  consequence 
of  this  focusing  asymmetry,  the  power  output  maximum 
occurs  at  a  frequency  slightly  below  gain  center.  Ex¬ 
perimental  verification  of  the  theory  is  described  in  the 
next  section. 

III.  Experiment 

The  gain  focusing  and  dispersion  focusing  have  been 
observed  experimentally  using  a  high  gain  3.51-/*  xenon 
laser.  In  a  discharge  laser  of  this  type  the  gain  maxi¬ 
mum  must  be  at  the  axis  of  the  discharge  with  the  gain 
falling  to  zero  at  the  tube  walls.  This  gain  profile 
makes  the  xenon  laser  appropriate  for  studying  focusing 

CII  CClSt 

®l)ot  size  *n  a  simple  plane  parallel  resonator 
filled  with  an  unsuturated  high  gain  medium  is  given 
by  Eq.  (20).  For  frequencies  near  line  center  Daw. son’s 
ntegral  F(x)  goes  to  zero  and  the  spot  size  is  simply 

w(x)  =  (8X*/]r0j)*-  (37) 


Experimental  investigations  of  this  resuii  have  been 
reported  previously.2  The  purpose  of  this  section  is  to 
consider  experimentally  the  more  subtle  frequency 
asymmetry  resulting  from  dispersion  focusing. 

1  he  apparatus  is  shown  in  I'ig.  3.  The  dc  discharge 
was  about  5.5  mm  in  diameter  and  1.1  min  length.  The 
right  mirror  was  highly  reflecting  and  could  be  trans¬ 
lated  longitudinally  by  means  of  a  motor  drive.  The 
cavity  length  was  1.29  m,  so  the  emptv  cavity  mode 
spacing  would  be  about  c/2L  =  116  MHz.  Synchro¬ 
nous  detection  was  used  to  improve  the  signal-to-noise 
ratio,  and  monoisotopic  xenon  was  used  to  prevent  un- 
iiccc-ssaf.v  asymmetries  in  tne  output.  The  xenon 
pressure  was  maintained  at  about  5  m  by  means  of  a 
liquid  nitrogen  trap.8 

A  U’wl  !  L"  J  Mu  ix***  otHVrft  for  (Lutafeiug 
cavity  length  (increasing  frequency)  is  shown  in  Fig.  4. 
The  laser  was  operated  very  near  threshold  and  the 
peaks  represent  successive  longitudinal  modes.  These 
peaks  are  to  be  compared  to  the  theoretical  curves 
shown  in  Fig.  2.  In  the  experimental  plot  there  is  a  dip 
in  the  output  power  on  the  high  frequency  side  of  the 
peak.  This  is  the  Lamb  dip  and  it  results  from  the 
interaction  of  the  left  and  right  traveling  beams  with 
atoms  which  have  zero  z-component  of  velocity.  Thus 
the  Lamb  dip  provides  a  convenient  indication  of  the 
frequency  x  =  0. 

Comparison  of  the  experimental  and  theoretical  plots 
show's  that  the  power  maximum  is  shifted  down  in  fre¬ 
quency  by  roughly  the  amount  predicted  by  the  dis¬ 
persion  focusing  theory.  Some  data  are  shown  in  Fig.  5. 
d  he  value  of  the  gain  as  a  function  of  discharge  current 
was  determined  by  introducing  known  losses  into  the 
cavity  and  reducing  the  current  until  threshold  was 
reached.  The  reflectivity  of  the  output  mirror  was  4% 
and  the  value  of  b  is  given  by  Eq.  (33).  For  the  analysis 
of  the  data  it  was  necessary  to  take  into  account  mode 
pulling,  because  with  a  dispersive  medium  of  this  sort 
the  rate  of  change  of  the  frequency  with  mirror  position 
may  be  greatly  reduced  near  line  center.9  Fora  rigor¬ 
ous  comparison  of  the  experimental  and  theoretical 
results  over  the  entire  spect!  ;m  it  would  have  been 
necessary  to  include  Ihe  nonlinear  mode  pulling  which 
occurs  near  the  wings  of  the  gain  line.  The  pressure 
was  low  enough  in  these  experiments  that  collision  ef¬ 
fects  are  believed  to  be  completely  unimportant.  The 
experiments  indicate  that  the  threshold  Lamb  dip  in 
xenon  has  a  width  of  about,  6  ±  1  MHz  and  a  depth 
somewhat  greatev  than  the  10%  enhanced  Lamb  dip 
reported  previously.10 

A  possible  cause  of  asymmetry  in  the  power  measure¬ 
ments  is  the  mass  motion  of  the  emitting  atoms.11 
Particularly  in  a  low  pressure  dc  discharge  one  might 
expect  a  drift  of  the  ions  toward  the  cathode  compen¬ 
sated  by  a  drift  of  neutral  atoms  toward  the  anode.  In 
a  high  gain  laser  the  Doppler  shifts  resulting  from  this 
mass  motion  would  result  in  an  asymmetry  of  the  out¬ 
put  power  spectrum.  However,  the  asymmetry  for 
light,  emerging  from  one  end  of  the  .  -  er  would  be  ex¬ 
pected  to  be  in  the  opposite  direction  to  the  asymmetry 
for  light  emerging  from  the  other  end.  To  check  this 


February  1972  /  Vol.  11,  No.  2  /  APPLIED  OPTICS  465 


Fig.  4.  Power  output  for  decreasing  cavity  length  with  a  dis¬ 
charge  current  of  18  mA. 


‘I 


Fig.  5.  Frequency  of  the  power  maximum  vs  the  threshold 
parameter  b.  The  solid  line  is  a  theoretical  result  obtained  from 
Eq.  (32)  and  the  dashed  line  is  a  plot  of  Eq.  (36).  The  circles 
are  experimental  values. 


possibility  a  resonator  was  constructedhaving  equally 
transmitting  mirrors  at  the  two  ends.  It  was  found 
that- the  power  spectrum  was  identical  at  the  two  ends 
of  the  laser,  indicating  that  for  the  conditions  of  our 
experiments  drift  of  the  atoms  is  unimportant.  The 
shift  resulting  from  the -known  abundances  of  impurity 
isotopes  is  estimated  to  be  small  compared  to  the  ob¬ 
served  asymmetry. 

IV.  Conclusion 

The  transverse  modes  of  a  laser  containing  a  high 
gain  medium  may  differ  significantly  from  the  modes  of 
a  similar  low  gain  laser  because  of  gain  and  dispersion 
focusing.  The  dispersion  effects  result  in  an  asymmetry 
of  the  output  power  spectrum  which  could  be  important 
in  any  Lamb  dip  measurements  in  lasers  with  moder¬ 
ately  high  gain.  Experiments  with  a  high  gain  3.51-p 
xenon,  laser  have  yielded  results  in  agreement  with  the 
theory.  While  only  plane  parallel  resonators  have 
been  considered  here,  the  results  may  readily  may  be 
extended  to  more  cor  'icated  laser  configurations. 

This  research  was  supported  by  the  Advanced  Re¬ 
search  Projects  Agency  through  the  Army  Research 
Office,  and  by  the  Air  Force  Office  of  Scientific  Re¬ 
search. 

References 

1.  H.  Kogelnik,  Appl.  Opt.  4,  1562  (1965). 

2.  L.  W.  Casperson  and  A.  Yariv,  Appl.  Phys.  Lett.  12,  355 
(1968). 

3.  A.  Szoke  and  A.  Javan,  Phys.  Rev.  145,  137  (1966). 

4.  B.  L.  Gyorffy,  M.  Borenstein,  and  W.  E.  Lamb,  Jr.,  Phys. 
Rev.  169, 340  (1968). 

5.  A.  Javan  and  P.  L.  Kelley,  IEEE  J.  Quant.  Electron.  QE-2, 
470  (1966). 

6.  L.  W.  Casperson,  Ph.D.  thesis,  California  Institute  of  Tech¬ 
nology  (1971),  Chap.  5. 

7.  D.  H.  Close,  Phys.  Rev.  153,  360  (1967),  Eq.  (44). 

8.  D.  R.  Armstrong,  IEEE  J.  Quant.  Electron.  QE-4, 968  (1968). 

9.  L.  W.  Casperson  and  A.  Yariv,  Appl.  Phys.  Lett.  17,  259 
(1970). 

10.  S.  C.  Wang,  R.  L.  Byer,  and, A.  E.  Siegman,  Appl.  Phys. 
Lett.  17, 120  (1970). 

11.  A.  D.  White,  Appl.  Phys.  Lett.  10,  24  (1967). 


SWfttJlitStesr- 


Ref.  2 

IEEB  JOUBNAL  OF  QUANTUM  ELECTRONICS,  VOL.  QS-8,  NO.  2,  FEBBUABT  1972 

Spectral  Narrowing  in  High-Gain  Lasers 
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Abstract — The  dependence  of  spectral  narrowing  in  lasers  on  the 
line-broadening  mechanism  is  investigated  including  the  effects  of 
saturation  and  distributed  loss.  It  is  found  that  in  the  unsaturated 
regime  the  narrowing  is  essentially  independent  of  the  resonance 
broadening  mechanism  and  the  narrowed  line  approaches  a  Gaus¬ 
sian.  The  onset  of  saturation  slows  or  reverses  the  narrowing 
process.  Experimental  results  have  been  obtained  using  a  3.51-m 
Xenon  laser. 


I.  Introduction 

SPECTRAL  narrowing  refers  to  the  fact  that  under 
some  circumstances  radiation  incident  on  a  laser 
amplifier  or  generated  within  an  amplifier  will 
emerge  with  a  spectrum  which  is  narrower  than  the  one  it 
started  with  [1],  In  high-gain  lasers  spectral  narrowing 
may  be  substantial  and  can  provide  a  highly  stable  and 
monochromatic  light  source.  The  applications  of  such 
optical  frequency  standards  in  metrology  are  well 
known  [2].  Comparison  of  the  input  and  output  spectra 
of  a  laser  amplifier  might  provide  a  sensitive  indirect 
measurement  of  the  amplifier  gain. 

The  purpose  of  this  paper  is  to  study  in  some  detail 
the  influence  of  an  amplifying  medium  on  a  spectral 
continuum  including  the  effects  of  saturation  and  dis¬ 
tributed  loss.  Emphasis  is  placed  on  the  important  prob¬ 
lem  of  supcrradiance,  but  narrowing  in  other  amplifiers 
is  also  considered  and  limiting  line  widths  are  determined. 

The  starting  point  for  this  investigation  is  a  general 
expression  for  the  incremental  gain  of  a  monochromatic 
signal  in  a  laser  amplifier.  From  the  derivation  of  Gor¬ 
don  et  al.  [3]  the  intensity  of  radiation  at  the  frequency 
fj  in  a  Doppler  broadened  medium  is  described  by 

d[i  _  Mi  T  _ exP  (-J.V) - =  dy 

*  [1  +  (2/  -  y.)J][ 1  +  1  +  (y  -  y]y\ 

(1) 

where  y(v)  =  2(v  —  v0)/Avk  is  a  normalized  frequency, 
k  is  a  pumping  constant,  s  is  a  saturation  parameter, 
and  the  natural  damping  ratio  t  =  (AvJ  AvD)  Vln  2 
measures  the  relative  importance  of  homogeneous  and 
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Doppler  broadening.  The  homogeneous  and  Doppler  line 
widths  are  given  respectively  by  Avk  and  Apd-  The  deriva¬ 
tion  can  be  generalized  in  a  straightforward  fashion  to 
the  interaction  with  an  optical  continuum  and  the  result  is 


dl(y i)  _  kl(y,) 

dz  r 


cxp(-cV) 


[l  +  (y  -  y\Y] 


[■+*/:. 


Hvn)  dy „  ~1 

l  +  (2/  -  yS  J 


dy  (2) 


where  I  (yi)  is  the  spectral  density  at  the  frequency  yi. 
This  result  neglects  any  coherent  interaction  between 
the  frequency  components  of  the  spectrum.  However,  it  is 
found  in  practice  that  even  with  monochromatic  fields 
coherence  effects  are  usually  unimportant  [4]  and  the 
neglect  of  coherence  makes  possible  a  number  of  useful 
analytic  solutions.  Also,  the  extreme  saturation  that  is 
necessary  for  coherence  effects  is  generally  found  to  im¬ 
pede  the  narrowing  process,  so  that  in  any  practical 
situation  strong  saturation  would  be  avoided. 


II.  Unsaturated  Amplifiers 

In  an  unsaturated  amplifier  the  growth  of  an  intensity 
continuum  is  governed  by 


~  =  —  aI(y'>  +  vg(yi)-  (3) 

The  unsaturated  incremental  gain  spectrum  g{y i)  is,  for 
simplicity,  assumed  to  be  independent  of  the  distance  z. 
The  second  term  on  the  right  side  of  (3)  represents  dis¬ 
tributed  losses.  The  last  term  is  the  spontaneous  “mis¬ 
sion,  which  has  the  same  frequency  dependence  g(y{) 
as  the  incremental  gain.  This  one-dimensional  model  is 
approximately  valid  in  narrow  bore  amplifiers  having 
a  length  much  greater  than  the  diameter.  The  coefficient 
rf  is  then  proportional  to  the  spontaneous  emission  rate 
and  to  a  geometrical  factor  that  depends  on  the  ampli¬ 
fier  dimensions.  A  more  rigorous  three-dimensional  treat¬ 
ment  would  have  to  account  for  the  spatial  distribution 
of  the  emitting  atoms  [5],  the  strong  focusing  effects  that 
are  common  in  high-gain  lasers  [6],  and  reflections  at 
the  boundaries  of  the  amplifying  medium.  Solving  for 
an  amplifier  of  length  z  yields 


I(y I,  -')  =  I(y  1,  0)  exp  [(0(2/,)  -  a)z] 

+  <exP  [&&*)  -  w 

Superradiance  will  be  considered  first.  In  a  super¬ 
radiant  source  there  is  no  input  and  if  losses  are  negli- 
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gible  (4)  simplifies  to  zero,  it  may  be  exparded  as 


I(y„  2)  =  >?(cxp  [giydz)  -  1).  (5) 

Defining  f(yt)  as  the  fraction  of  the  line  center  intensity 
at  the  frequency  yt  yields 


//,.  -1  =  KVu  2)  _  exp  [g(y,)z]  -  1 
ny"Z>  7(0,  z)  exp  [<?(0)z]  —  1 


(6) 


Most  lasers  can  be  classed  as  either  homogeneously  or 
inhomogeneously  broadened  and  only  these  cases  will  be 
considered  here. 

For  homogeneous  broadening  (e  »  1)  and  no  satura¬ 
tion  the  gain  from  (2)  is  simply  the  Iorentzian 


Q(lJl)  bom  /—  |  ,  2*  (7) 

V  IT  t  1  +  2/l 

The  spectral  width,  defined  as  the  separation  between 
the  two  frequencies  at  which  the  spectral  intensity  is 
down  to  \  its  peak  value,  is  obtained  by  combining  (7) 
with  (6)  for  /  =  ^.  The  result  is 

Aph„ra  =  — - - ~-7=~, - 1.  (8) 

*  Vt  c  In  j (exp  {kz/Vrt  ]  +  1) 

For  short  distances  ( kz /  \/ir<  «  1)  the  line  width  given  by 
(8)  is  just  the  homogeneous  line  width  Ac*.  For  long 
distances  (kz/  >>  1)  (8)  simplifies  to 

Aphom  A.^  In  2  =  ^lg(Q)b~z'  (9) 

Thus  the  narrowing  effect  becomes  important  when  the 
product  {/(O)hOTi  2  becomes  comparable  to  unity.  The  ap¬ 
proximate  result  given  by  (9)  becomes  valid  after  the 
width  is  narrowed  to  about  one  half  of  its  initial  value. 

For  an  unsaturated  inhomogeneously  broadened  am¬ 
plifier  (t  «  1)  and  for  radiation  not  too  far  in  the  wings 
of  the  line  (yt  <K  1/2  e2)  the  gain  from  (2)  is  the  Gaus¬ 
sian 


siyO  —  ffo  —  (hy/  +  +  <7»J/i4  +  •  •  •  (13) 

where  g0  and  g2  are  positive.  The  spectrum  of  the  super¬ 
radiance  is  given  by  (6)  as 

Kvi) 


=  exp  f(<7„  -  gd/j*  4-  g3?/|3  +  •  •  • )z ]  -  1 
exp  [go2]  -  1 


exp  \gnz)  -  1 

where  A„  = 

At  large  distances  ( g0z  »  1)  the 
cancel /leaving 


exp  [(  )']  •  ••  -  1 
(14) 

zero-order  terms 


f(y,)  ~  exp  [-(^)  ]  exp  [(■&)  ]  exp  [(  )4]  -  •  ••  .  (15) 

Also,  at  very  large  distances  one  finds  that  A„  «  A„,j 
so  that  the  Gaussian  factor  is  much  narrower  than  the 
others.  Consequently,  all  of  the  factors  but  the  first  may 
be  replaced  with  their  value  at  line  center.  Therefore, 

/(»,)=,  «p  [-(*)’]  ,  (16) 

This  is  in  agreement  with  the  previous  results  for  homo¬ 
geneous  and  inhomogeneous  broadening  as  may  be  veri¬ 
fied  by  expanding  g(yt)  of  (7)  and  (10)  in  power  series 
in  yi  to  obtain  g2.  More  generally  if  the  incremental  gain 
has  several  maxima,  the  emission  will  eventually  resolve 
itself  into  narrowing  Gaussian  lines  centered  on  the  gain 
maxima.  This  resolving  effect  has  been  observed  by 
Parks  et  al.  [7]- 

So  far  only  superradiance  has  been  considered.  If  an 
amplifier  has  an  input  I(yi,  0)  and  negligible  spontane¬ 
ous  emission,  (6)  for  the  output  light  spectrum  is  re¬ 
placed  by 


g(>J ()ui,on,  -  k  exp  (~t\ /,2).  (10) 

Combining  this  with  (6)  for  /  =  $  yields 

.  ,  /in  kz  -  In  In  £[exp  (kz)  +  l]  /t  ,N 

ACjohom  =  A  VDyJ - \/l2  - -  '  (U) 

For  short  distance  (kz  1)  the  width  of  the  emission 
is  equal  to  the  Doppler  width  Apc.  For  long  distances 
(kz  »  1)  one  finds 


Acin 


Ar0  7=-  = 


Ar0 


V  kz  V  g(  0)  inhom^ 


(12) 


These  results  suggest  that  the  narrowing  proceeds  in 
about  the  same  fashion  independent  of  the  line  broad¬ 
ening  mechanisms.  To  verify  this  conclusion,  one  may 
consider  a  completely  general  incremental  gain  function 
g(iji)  with  a  maximum  at  the  frequency  y0  =  0.  Assum¬ 
ing  that  g(yi)  is  differentiable  in  the  neighborhood  of 


tv,.  _  HVi,  0)  exp  \g(ih)z] 
KV"  ]  7(0,0)  exp  [gWzj ' 


At  large  distances  this  becomes 


F(l l„z) 


7 (y,,  0) 
7(0,  0) 


exp 


[-©’]  • 


Ai  = 


V  g2z 


(17) 


(18) 


If  the  input  is  reasonably  smooth  the  output  spectrum 
is  simply 


F(>Ji ,  2)  exp  £  -  J  •  (19) 

If  the  input  is  a  narrow  Gaussian  of  width  Aln,  the  out¬ 
put  will  be  a  Gaussian  of  width  A™,  such  that 


(£)’  -  (iJ + (£)’■  <20) 


The  concept  of  gain  narrowing  is  sometimes  useful 
instead  of  spectral  narrowing.  The  two  are  obviously 
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clo-elv  related.  If  one  that  tlic  input  spectrum 

in  1 17)  i:-  white.  t lion  the  intensity  factors  cancel  and 
Fi  iji)  is  a  general  expression  for  the  gain  spectrum  of 
an  amplifier  normalized  to  unity  on  line  center.  Expres¬ 
sions  for  the  gain  line  width,  for  example,  are  then  ob¬ 
tained  by  setting  /•'  equal  to  one  half  and  solving  (17) 
for  ?/n/2-  For  an  inliomogcneotisiy  broadened  amplifier 
the  results  may  be  made  to  conform  with  those  of 
Hotz  [8]. 

In  summary,  one  may  conclude  that  in  an  unsaturated 
amplifier  the  width  of  the  amplified  spectrum  decreases 
with  distance.  The  shape  of  the  spectrum  approaches 
a  Gaussian  as  the  effective  part  of  the  gain  approaches 
a  quadratic.  As  an  example,  in  a  helium-xenon  discharge 
the  gain  may  be  400  dB/m  [o|  or  nearly  k  —  100  so  that 
the  spectrum  of  a  superradiant  source  1  m  long  would 
be  narrowed  according  to  (12)  by  a  faetor  of  ten  pro¬ 
vided  saturation  did  not  occur.  The  Doppler  width  of 
xenon  at  room  temperature  is  about  Arn  — 1  100  MHz, 
so  the  narrowed  radiation  would  have  a  width  of  about 
10  MHz.  Since  the  frequency  of  the  3.ol-/e  transition 
is  about  10’  MHz,  the  light  would  be  monochromatic  to 
a  part  in  107.  Collimation  in  a  long  high-gain  gas  laser 
is  taken  eare  of  by  the  gain  profile  of  the  medium  it¬ 
self  [G]  and  there  should  not  be  much  difficulty  in  con¬ 
structing  a  gas  discharge  amplifier  of  arbitrary  length. 
However,  care  is  necessary  to  prevent  saturation  as  will 
be  shown  in  the  next  scetion. 

III.  Saturation  Effects 

The  intensity  of  amplified  spontaneous  emission  can 
easily  reaeh  the  saturation  level.  When  this  happens  the 
behavior  of  the  spectrum  becomes  considerably  more 
complicated  and  for  a  completely  general  treatment 
computer  solutions  are  required.  We  present  here  ana¬ 
lytic  solutions  for  certain  important  limiting  situations 
and  also  some  general  computer  results. 

For  a  homogeneously  broadened  amplifier  all  of  the 
factors  but  the  Gaussian  may  be  removed  from  the 
integral  in  (2).  The  result  of  the  integration  is 

dl(y.)  =  kl  (?/,)  I 

dz  vV  fh(z)  1  -f*  Vi 

where 

1  /  \  ,  I  ,  r  I(jfn)  d]l„ 

and  distributed  losses  are  ignored.  Thus  the  gain  profile 
remains  Lorentzian  even  with  saturation,  although  its 
amplitude  decreases.  The  quadratic  term  in  the  frequency 
expansion  of  the  gain  decreases  in  magnitude  and  hence 
a  narrow  Gaussian  beam  will  continue  to  narrow,  but 
at  a  redueed  rate  as  saturation  becomes  important.  Solv¬ 
ing  as  before  yields  for  large  distances 

A^bom  =  Vt  €  In  2  j k  J  '  (23) 


(21) 

(22) 


To  proceed,  expressions  for  h(z)  must  be  obtained. 

For  an  intensity  spectrum  that  is  narrow  compared 
to  Ai'„  (22)  becomes 

h(z)  ~  1  +  s  [  I(y„)  dyn  =  1  +  si,  (24) 

where  I,  is  the  total  intensity,  Such  a  speetrum  would 
be  obtained,  for  example,  at  large  distances  ( g  (0)  flr.m2 
3>  1)  in  a  superradiant  amplifier.  To  evaluate  h{z), 
it  is  first  necessary  to  find  expressions  for  I,.  In  this  ap¬ 
proximation  (21)  may  be  integrated  over  frequency 
yielding 


dl,  =  U, _ 

dz  v^t(l  +  si,) 


(25) 


Thus  in  an  unsaturated  amplifier  the  intensity  grows 
according  to 


1,  =  I,0  exp  [b/Vre]  (26) 

while  in  a  highly  saturated  amplifier  the  intensity  is 
governed  by 


/,  =/,  o  +  -S—  (27) 

V7T  tS 

The  intensity  in  the  saturated  amplifier  eventually 
reaches  the  point  at  which  the  loss  term  —  alt,  which 
was  neglected  initially,  is  comparable  to  the  saturated 
gain  term.  If  losses  are  included,  the  intensity  is  gov¬ 
erned  by 


dh  =  kl, 
dz  -\Ztt  c(l  +  si,) 


—  al,  =  0, 


which  for  si,  »  1  has  the  steady-state  solution 


(28) 


/,  =  -  (29) 

V  tt  eas 


These  results  may  be  collected 


lossless 

unsaturated 

regime 

lossless 

saturated 

regime 

loss-limited 

saturated  regime.  (30) 


They  are  similar  to  Rigrod’s  solutions  [9]  for  monochro¬ 
matic  radiation  at  gain  center. 

Using  (24)  and  (30),  (23)  may  be  written  for  large 
distances  as 
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lossless 
unsaturated 
regime  « 

lossless 
saturated 
regime 

loss-limited 

saturated  regime.  (31) 

These  results  are  valid  provided  that  the  spectrum  is 
much  narrower  than  the  homogeneous  line  width  before 
saturation  sets  in.  As  the  gain  is  pulled  down  by  satura¬ 
tion,  the  narrowing  rate  is  slowed.  When  losses  become 
important,  the  gain  curve  is  elamped  and  the  narrowing 
speeds  up  again.  In  a  homogeneously  broadened  amplifier 
neither  saturation  nor  losses  stop  the  narrowing  process. 
In  Fig.  1  are  some  numerical  solutions  for  the  line  width 
of  the  emission  from  a  homogeneously  broadened  super¬ 
radiant  laser  amplifier  as  a  function  of  the  normalized 
distance  Zhom  =  kz/(  y/vt  In  2).  For  simplicity  losses 
are  assumed  to  be  negligible.  The  parameter  in  these 
plots  is  the  product  S  =  sij. 

For  an  inhomogeneously  broadened  amplifier  (t  «  1) 
it  will  be  assumed  that  the  intensity  is  nearly  uniform 
over  a  natural  line  width.  Then  the  intensity  spectrum 
may  be  removed  from  the  denominator  integral  in  (2) 
and  the  result  simplifies  to 

i  1  dl(y ,)  k  exp  [-tWl  /32) 

I(yi)  dz  1  +  irsl(iji) 

The  speetral  density  ,r(yi)  is  found  in  the  various  regions 
from  (32)  in  a  mannei  essentially  identical  to  the  homo¬ 
geneous  ease.  The  results  are 


Arh. 
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V4 


In  2 


In  2 


'In 


kz 

VV 
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In  2 
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Hy  i, *) 


lossless 

I(y,,  0)  exp  [kz  exp  [— «2y42]l 
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regime 

lossless 

I(Vi,  0)  +  —  exp  [-«Vl 

■JTO 

saturated 

regime 

loss-limited 

~~  exp  [-«y  ] 

saturated 

iraS 

regime 

and  consequently  the  speetral 


(33) 

width  for  large  distances  is 


Fig.  2.  Superradiant  narrowing  in  a  Doppler  broadened  amplifier 
for  various  values  of  8  =  sij. 

It  is  evident  that  the  cffeets  of  saturation  on  narrow¬ 
ing  in  an  inhomogeneously  broadened  amplifier  are  sig¬ 
nificantly  different  from  the  effects  in  a  homogeneously 
broadened  amplifier.  In  the  inhomogeneous  ease  the  onset 
of  saturation  reverses  the  narrowing  process  and  restores 
the  radiation  to  its  Doppler  line  shape.  This  oecurs  be¬ 
cause  the  center  of  the  line  saturates  first,  while  the  wings 
eontinuc  to  grow  exponentially.  This  analysis  is  valid 
provided  that  the  intensity  speetrum  remains  broad  com¬ 
pared  to  the  homogeneous  line  width.  If  the  spectrum  be¬ 
comes  narrow  compared  to  Ar»,  the  rebroadening  would 
be  expected  to  occur  more  slowly.  Some  numerical  re¬ 
sults  for  superradiant  narrowing  in  an  inhomogeneously 
broadened  laser  are  shown  in  Fig.  2. 

The  minimum  line  width  for  a  simple  inhomogeneous 
superradiant  source  may  be  readily  ealeulated.  From 
(32)  it  is  evident  that  saturation  becomes  important 
when  I(y0)  =  1/n-s.  Then,  using  (5),  one  finds  that 
saturation  oeeurs  at  a  distance  zeBt  given  approximately 
by 


I 

y/kz 

Ap  j  nhom  _ 

Apd 

1 


lossless 

unsaturated 

regime 

lossless 

saturated 

regime 

loss-limited 
saturated  regime. 


~  —  In  tStj. 


(35) 


Use  of  this  expression  in  (12)  yields  the  minimum  line 
width 


(34) 


Ap 


i  ohom  min 


_ Ap d__ 

y/  —  111  irStj 


(30) 
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.Thus,  if,  for  example,  ,  ^  5  X  llh«  W/m*  and  s  ^  6  X 
10~»  m’/w  in  a  xenon  laser,  then  the  Doppler  line  can  be 
narrowed  by  at  most  a  factor  of  about  5.5  before  satura¬ 
tion  becomes  important  at  a  distance  (for  k  ~  15)  0f 
about  2.0  m.  These  values  for  v  anc.  s  are  judged  to  be  rea¬ 
sonably  valid  for  a  simple  xenon  laser  used  in  some  of  out 
experiments  based  on  the  data  of  Clark  [10].  This  nar¬ 
rowing  is  not  too  impressive,  and  deereasing  v  by  orders 
of  magnitude  does  not  help  much,  sinee  the  dependence 
on  v  involves  a  logarithm  and  a  square  root. 

A  possible  scheme  for  reducing  the  ultimate  line  width 
is  to  place  attenuators  between  sections  of  the  amplifying 
medium  [11].  These  would  cut  down  the  intensity  to 
prevent  saturation  without  affecting  the  narrowing  proe- 

ess,  iEven  in  such  a  system,  however,  the  spectral  line 
width  could  never  approach  zero  because  there  is  always 
broadband  background  noise  being  added  to  the  beam 
by  spontaneous  emission.  The  result  of  the  background 
is  that  the  spectrum  must  eventually  approach  a  narrow 
limiting  line  shape. 

The  narrowest  possible  line  would  be  obtained  in  a 
long  amplifier  with  distributed  losses  which  are  just  suf- 
ficient  to  keep  the  line  center  intensity  somewhat  below 
the  saturation  intensity  1  As.  To  get  an  estimate  of  this 
limiting  line  shape  one  can  write  (3)  for  steady  state 
with  /(yo)  =  1  /-ns 


Using  the  approximate  numbers  given  previously  for  s 
and  v,  one  finds  that  the  Doppler  line  would  be  nar¬ 
rowed  by  a  factor  of  about  4  x  IQ-7.  A  Doppler  width 
of  aVd  ~  10*  Hz  could  yield  an  intensity  spectrum  of 
a  out  40-Hz  width.  If  the  oscillation  frequency  were 
about  10“  H"  as  in  xenon,  the  output  could  be  used  as 
an  absolute  frequency  standard  with  a  stability  of  about 
four  parts  m  lO13.  Similar  calculations  can  be  carried  out 
for  the  limiting  line  shape  in  a  laser  incorporating  dis¬ 
crete  rather  than  continuous  losses. 

The  preceding  discussion  suggests  that  superradiant 
lasers  could  be  useful  as  extremely  stable  frequency 
standards.  Some  practical  limitations  on  such  a  system 
should  be  emphasized.  The  intensity  only  approaches  its 
limiting  form  at  a  rate  given  by  (12).  Thus,  to  obtain  a 
me  width  of  40  Hz  for  a  gain  constant  of  k  ~  100  m-1 
the  overall  length  of  the  laser  would  have  to  be  greater 
than  10s  m.  However,  higher  gain  media  are  available 
and  for  some  applications  superradiant  lasers  should  be 
useful  as  absolute  frequency  standards. 


n  =  £°  _  SL  i  „ 

°  xs  rs+  v9°- 


(37) 


Thus,  the  appropriate  value  for  the  loss  constant  «  is 

a  =  0o(l  +  tst,).  (38) 

Using  this  result,  (3)  away  from  line  center  can  be 
written  at  steady  state  as 

0  =  a(y,)i(y,)  -  00(1  +  *sv)I(y,)  +  vg(y,)  (39) 
with  the  solution 


Hyd  = 


vg(y<) 


g»( 1  +  —  g(y,) 


1  +  rSr)  — 


g(yi) ' 


(40) 


go 


Keeping  the  second-order  term  in  the  power-series  ex¬ 
pansion  of  g< ?/,)  leads  finally  to  the  intensity  spectrum 


7(i/,)  - - 

1  +  ( gt/rsijg0)yi * 


(41) 


Therefore,  the  narrowest  possible  line  is  a  Lorentzian  of 
width 


i  n  — 


fggggn 

g* 


(42) 

If  the  gain  profile  is  the  Gaussian  given  by  (10)  then 
02  —  0ov  and  the  line  width  is  simply 


*  Vx Si]  =  . 


(43) 


IV.  Experiment 

We  describe  here  an  experiment  that  has  been  per¬ 
formed  in  an  effort  to  verify  some  of  the  conclusions  of 
c  previous  sections  regarding  spectral  nar-owing.  The 
apparatus  used  for  this  study  consisted  of  a  dc  xenon 
discharge  having  an  active  region  of  5.5  mm  diameter 
and  U-m  length.  The  xenon  pressure  was  maintained 
at  aoout  5  /*  by  means  of  a  liquid  nitrogen  trap  [12].' 

source8861,  ^  °Perated  85  a  Sing,e  mirror  superradiant 

To  measure  the  width  of  the  chopped  superradiant 
output,  we  allowed  it  to  impinge  on  an  InAs  junction 
detector  The  resulting  current  was  then  fed  into  a  con¬ 
ventional  RF  spectrum  analyzer  followed  by  a  lock-in 
amplifier.  It  can  be  shown  that  under  these  conditions  a 
uaussian  optical  intensity  spectrum  7(f)  will  -ive  rise 
to  a  low  frequency  spectrum  [13]. 

«  l  /({)/({  +  v)  (44) 

tah«ha.t,°7«rian  and  Wh°Se  "idtt  ^  ^  by  * 

nflntyPiCaUPeCtrUm  13  shown  in  Fjg-  3  for  a  current 
l  P13  CUrrent  corresPouds  to  a  gain  of  about 

?  ~  U  m.  *  The  spectrum  is  approximately  Gaussian 
in  shape  with  a  width  of  about  20  MHz.  Therefore,  the 
width  of  the  original  intensity  spectrum  is  about  14  MHz 

is  2  2e  °f  the  la3er  amPlif>^g  medium 

2.2  m.  Using  (12),  the  intensity  spectrum  should  be 

nan-owed  from  the  Doppler  width  by  the  factor  \/fz  = 

lOO  MHz^tr  5i4‘  -f  the.D°PPler  ™dth  were  AvD  = 
iq°^  Vm  thlen,  he  mtensity  spectral  width  should  be 
19.5  MHz  which  is  in  satisfactory  agreement  with  the 

Z71  r  T  °f  H  ^HZ>  The  discrepancy,  if  significant, 

oould  be  due  to  a  small  amount  of  nonresonant  feedback 
scattered  back  into  the  laser.  Care  was  also  necessary  to 
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Fig.  3.  Superradianl  output  spectrum  for  nn  SO-mA 
discharge  current  with  5-MIIz/div  dispersion. 


prevent  resonant  feedbaelc,  which  eould  result,  for  ex¬ 
ample,  from  reflections  off  the  deteetor.  Resonant  feed¬ 
back  led  to  longitudinal  modes,  whieh  could  be  detected 
as  periodic  intensity  fluctuations  as  the  single  mirror  was 
seanned  longitudinally. 

Some  experimental  data  are  eolleeted  in  Fig.  4.  Accord¬ 
ing  to  (12)  the  superradiant  line  width  should  be  in¬ 
versely  proportional  to  the  square  root  of  the  gain  con¬ 
stant,  which  is  in  satisfactory  agreement  with  the  data. 
Line  width  data  have  not  been  obtained  for  lower  values 
of  gain  beeause  of  the  limited  sensitivity  of  our  detection 
system.  At  higher  levels  of  gain  the  laser  amplifier  in¬ 
evitably  began  to  oseillate  so  that  again  meaningful  data 
could  not  be  obtained.  These  are  only  preliminary  results 
and  more  refined  experiments  are  in  progress. 

V.  Conclusion 

It  has  been  shown  in  this  paper  that  a  gain  profile 
whieh  is  quadratic  in  frequency  near  its  maximum  ean 
support  a  narrowing  Gaussian  radiation  speetrum.  The 


spectral  width  varies  inversely  as  the  square  root  of  dis¬ 
tance  at  long  distances.  In  a  homogeneously  broadened 
amplifier  saturation  slows  the  narrowing  process,  while 
in  an  inhomogoncously  broadened  amplifier  saturation 
restores  the  line  to  its  original  inhomogeneous  line  shape. 
Narrowed  amplifiers  are  useful  in  spectroscopy,  gain 
measurements,  and  anywhere  a  stable  absolute  frequency 
standard  is  needed. 

Analogous  results  are  obtained  for  laser  oscillators  if 
the  loordinate  z  is  replaced  by  the  time  coordinate  ct/no. 
One  finds  that  the  spectrum  of  an  abruptly  started  laser 
oscillator  should  be  a  narrowing  Gaussian  that  ap¬ 
proaches  ultimately  a  narrow  limiting  Lerentzian  line 
because  of  spontaneous  emission.  The  nat  rowing  in  this 
ease  results  from  the  Fabry-Perot  transmission  reso¬ 
nances  of  the  oseillator.  Nonresonant  osc'llators  operated 
near  threshold  may  also  be  useful  as  freqi  eney  standards. 
Nonr>sonance  ean  be  obtained  by  making  fn,»  cavity  long 
so  that  the  modes  are  closely  spaced  or  by  replacing  one 
of  the  eavity  mirrors  with  a  scatterer  [14]. 
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The  nonlinear  coefficients  of  Ki)l»  Hi) a  ami  I  iin 


We  wish  to  report  the  measurement  of  the  non¬ 
linear  coefficients  of  RbH2P04(RDP),  RbH2As04 
(RDA)  and  LiI03  relative  to  r/36  of  KJl2P04(KDP) 
[1|T.  Previous  measurements  [2|  of  1/35  of  RDP 
and  RDA  used  pulsed  laser  sources  and  thus  are 
subject  to  a  higher  degree  of  uncertainty  than 
measurements  employing  cw  sources.  The  non¬ 
linear  coefficient  of  LUO3  has  previously  been 
measured  to  good  accuracy  [3|  and  was  included 
in  our  measurements  as  a  check  on  the  accuracy 
of  the  experimental  technique. 

Our  experimental  setup  employed  a  cw  ruby 
as  the  fundamental  light  source.  A  schematic  of 
ihe  experimental  arrangement  is  shown  in  fig.  1, 
The  laser  is  similar  to  that  described  by  Evtuhov 
and  Neeland  [4|  and  operated  at  an  output  power 
level  of  100-200  mW.  The  output  power  was  con¬ 
tinuously  monitored  and  an  iris  inside  the  cavity 
assured  TEMqq  operation.  The  second  harmonic 
power  at  3412  A  was  detected  using  a  vacuum 
photodiode  and  standard  lock-in  amplifier  tech¬ 
niques. 

Because  of  tne  rapidly  diverging  output  of  our 
laserft  and  in  order  to  achieve  easily  observable 
UV  outputs,  it  was  convenient  to  focus  tightly 
into  the  nonlinear  crystals.  Two  different  focal 
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t  d-jcfKDP)  (1.0  ±0.4)  10-9  esu. 

I  he  output  confoeal  parameter  of  the  laser  was 
measured  (5)  to  be  h  2,3  cm  where  b  -  27(±)2/A  am| 
w0  is  the  minimum  gaussian  beam  radius.  0 
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I1  ig.  1 .  F.xperimental  arrangement  lor  measuring  rehi- 
tive  second-harmonic  powers  and  focused  beam  pro- 
liles.  Laser:  e u  rubj  :  BS:  beam  splitter:  f] :  focusing 
lens,  la:  collecting  lens:  PI  and  P2:  crossed  polarizers: 

I- :  red-cut.  1  V-pass  niters:  C :  eho|)per:  S:  73  p  slit: 

1)1:  phototube  with  S-4  surface:  1)2:  phototube  with  S-i 
surface. 

length  lenses  were  used  and  the  transverse  pro¬ 
file  of  the  focused  beam  scanned  with  detector- 
si  it  combination  to  obtain  the  confoeal  parameter 
in  atrj.  The  measured  confoeal  parameters  were 
then  used  with  the  measured  relative  UV  powers 
and  the  focused-beani  theory  of  Boyd  and  Klein- 
111311  I6  I  to  obtain  the  nonlinear  coefficients  rela¬ 
tive  to  r/36(KDP). 

For  KDP,  RDP  and  RDA.  all  losses  are  as¬ 
sumed  to  be  negligible  at  both  6943  A  and  3472  A 
The  loss  at  6943  A  in  LiI03  is  also  neglected, 
but  the  crystal  used  in  our  measurements  does 
have  a  finite  absorption  at  3472  A.  On  the  basis 
of  measurements  on  a  number  of  crystals  [7],  we 
have  used  o  ~  2  enr1  for  the  extraordinary  second- 
harmonic  wave  in  our  LiIQ3  crystal. 


f  The  confoeal  parameter  of  a  propagating  gaussian 
beam  can  be  determined  by  measuring  the  spot 
radius  at  tv  o  points  along 'the  direction  of  propaga¬ 
tion  and  applying  the  formulas  of  Pearson  et  al.  [ 5 j . 
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Table  1 

Experimental  values  and  relative  nonlinear  coefficients 


Crystal 

Phase 

matching- 

angle 

Crystal 
length  (cm) 

ci 

d/d' j(j(KDP) 

previous  |ref.J  now  value 

Kl)l> 

51° 

.’.3 

(/36 

1.0 

1.0 

KDP 

07° 

1.8 

(/30 

1.04  ±  0.15  [2] 

0.92  ±  0.1 

It  DA 

90° 

1.9 

30 

0.04  1  0.15  (2 J 

1.04  i  0.1 

UK  >3 

52° 

0.4 

'Da 

11  t  1.3  [3J 

11.2  ±  1.2 

The  results  of  our  measurements  are  shown 
in  table  1  along  with  previously  reported  values 
for  the  relative  nonlinear  coefficients.  The 
values  of  d  dgg(KDP)  for  RDP  and  LiIO,3  agree 
with  previous  measurements  to  within  the  experi¬ 
mental  error,  but  the  value  for  RDA  is  70°?  larger. 
All  crystals  were  angle-tuned  except  RDA  which 
was  temperature-tuned.  For  our  RDA  crystal, 
manufactured  by  Quantum  Technology,  phase¬ 
matching  at  90°  to  the  optic  axis  occurred  at  a 
crystal  temperature  of  97.4°C  and  the  second- 
harmonic  versus  temperature  scan  had  a  full- 
width  at  half-maximum  of  approximately  1.8°C. 

The  authors  would  like  to  express  their  in¬ 
debtedness  to  Dr.  Viktor  Evtuhov  for  the  loan  of 


the  KDP,  RDP,  and  LUO3  crystals  and  for  his 
h'.-ipful  comments  and  suggestions. 
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* 

A  proposal  for  a  new  method  of  phase  matching  in  nonlinear  optical  Interactions  Is  made.  A 
periodic  perturbation  of  the  surface  of  a  thin-film  waveguide  generates  space  harmonics 
with  new  propagation  constants  which  can  be  phase  matched.  An  analysis  of  this  proposal 
shows  it  to  be  particularly  interesting  for  a  class  of  thin-film  nonlinear  devices  using  the 
cubic  optically  isotropic  semiconductors  (such  as  GaAs,  GaP,  etc.)  which  possess  high  non¬ 
linear  optical  coefficients  but  are  not  phase  matchable  by  the  conventional  birefringent 
techniques 


Optically  isotropic  materials  have,  to  date,  not  been 
used  in  nonlinear  optical  applications  because  the  con¬ 
ventional  technique  of  birefringent  phase  matching  can¬ 
not  be  applied.  This  has,  so  far,  precluded  the  use  of 
many  highly  nor'inear  materials,  such  as  GaAs,  from 
practical  utilization  in  second-harmonic  generation, 
parametric  oscillation,  and  frequency  upconversion. 

One  approach  to  phase  matching  in  optically  isotropic 
materials  involves  the  use  of  dimensional  dispersion  in 
thin-film  waveguides.1  Another  suggestion2'3  utilizes 
periodically  laminated  structures.  The  realization  of 
this  last  approach  involves  a  fractional  wavelength  con¬ 
trol  of  the  lamination  period  and  has  not  yet  been 
demonstrated. 

As  an  outgrowth  of  our  experiments  on  light  guiding  in 
epitaxial  GaAs  thin  films,4  we  have  considered  a  new 
spatial  modulation  approach  to  phase  matching,  which 
involves  a  periodic  corrugation  of  one  (or  more)  of  the 
boundaries  of  the  tin. -film  waveguide,  as  shown  in 
Fig.  1. 

A  more  careful  investigation  of  this  idea  leads  to  the 
conclusion  that  the  expected  magnitude  of  the  effect, 
using  presently  available  materials  and  techniques,  is 
such  as  to  encourage  serious  efforts  to  implement  it. 
The  main  results  of  the  analysis  are  presented  in  what 
follows . 

Consider  a  single  propagating  mode,  say  m,  in  a  thin 
(uncorrugated)  waveguide  with  a  principal  transverse 
field  component, 

£"(*,  z ,  1)  =  C“  exp[i(wf  -  0  "z)]«?“(x) ,  (1 ) 

in  the  presence  of  a  sinusoidal  boundary  perturbation 
whose  period  A  is  large  enough  so  as  not  to  couple  into 
the  continuum  (“leaky”)  modes5;  the  wave  has  the  Flo- 
quet  form 

K(x,  M  =  Cm(z)exp[iM  -  0“z)]«?“(.v) 

■o 

=  exp (i(w/-0"2)l«?"(.v)  JCA^exp[-  in(2n/\)z], 

(2) 

where  Cm{z)  is  periodic  in  A.  The  mode  consists  of  an 
infinite  number  of  space  harmonics,  each  with  its  phase 
constant 

0“B  =  0“+w2ff/A,  n  =  ±l,2, . . .  .  (3) 

Phase-matched  interactions  are  thus  no  longer  limited 
to  the  principal  value  of  0  but  can  involve  the  space 
harmonics.  As  an  example,  second-harmonic  genera¬ 
tion  can  be  achieved  by  matching  the  fundamental  (w=0) 


space  harmonic  at  ui  to  the  first  (n  =  ±l)  space  harmonic 
at  2w,  so  that 

0§“  =20“  ±  (2ir/A).  (4) 

The  penalty  for  using  the  space  harmonics  is  that  the 
conversion  efficiency  is  reduced  relative  to  the  phase- 
matched  interaction  in  the  bulk  by  a  factor  which  in  the 
example  just  quoted  is  approximately  equal  to 
lA^  I2.  A  meaningful  evaluation  of  the  feasibility  of 
phase  matching  by  periodic  surface  perturbation  pre¬ 
requires  a  solution  for  the  amplitudes  of  the  space 
harmonics  Amrl. 

A  solution  of  Maxwell’s  equations  for  a  TE  mode,  to  be 
described  elsewhere,  gives 

^2-2 U«(*)l i*1?{x,z)dx,  (5) 


where  A n2{x,z)  is  the  deviation,  due  to  surface  corruga¬ 
tion,  of  the  actual  index  (squared)  of  refraction  from 
that  of  a  planar  boundary  film.  k0~2n/\  and  &m(x )  is 
normalized  according  to  /.”  I  <?m(x)l2  <fx=  1 .  When  a,  the 
corrugation  amplitude,  is  small  compared  to  the  thick¬ 
ness  t,  we  can  replace  <?m(x)  in  Eq.  (5)  by  <?„(()),  and 
assuming  d2Cm/dz2  <<  PmdCjdz ,  we  obtain 


Cm(«)  =  Cm(0)exp{i 


fe£lg,(0)l2(»g-l)gA 

4tt0„ 


The  corrugation  thus  causes  a  (spatial)  phase  modulation 
of  the  mode,  so  that,  using  Eq.  (2)  and  a  Bessel-func- 


n 


l 


**  -1- 


n3  substrate 


FIG.  1.  Model  of  a  thin-film  waveguide  with  a  periodic 
perturbation  of  one  interface. 
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(13) 


tion  expansion  of  Eq.  (6),  we  get  the  following  solution  Pttu(l)  _  2u>zdhtlz  /Po\3/2 

for  the  mode:  •  P“  “  (»")W"  \«0  /  m  ' 


Ek{x, z,t)  =  CJO) exp [i(ut -  <f>)  ]<?"(*) 

m 

X  lDm  exp{-  j'[/3"  -  w(2it/A)]z}, 


where 

,r-fcgl^(0)l2(«I-l)aA  (8) 

■  4it0m 

Equation  (7)  gives  explicitly  the  relative  amplitudes  of 
the  space  harmonics  generated  by  the  surface  corruga¬ 
tion.  This  result  can  now  be  used  in  analyzing  nonlinear 
interactions  in  a  thin  film.  To  be  specific  we  consider 
an  example  of  second-harmonic  generation  in  which  both 
the  injut  (w)  and  output  (2u>)  are  in  the  same,  say  m=  0, 
waveguide  mode.  We  thus  have 


£"(*,  z  ,t)  =  2./  An  exp[i  (wt  -  0")2  l«?"(x) , 
E*u(x,  z ,  t)  =  Tj  Bm  exp[i'(2 w  t  -  02m")z  ]<?;;“(*) , 


0) 


where  0"  =  0"-m2it/A,  and  An  and  Bm  are  defined  by  Eq. 
(7).  The  second-harmonic  polarization  generated  by 
g0(x,z,t)  is  taken  as 


/>*"(*,  z ,  f)  =  dE  E  AnA,  exp{f[2wf  -  (0"  +  0?)z  ]}[<?"(*)  I2, 

(10) 


where  d  is  the  appropriate  bulk  nonlinear  tensor  ele¬ 
ment.  The  rate  of  growth  of  the  average  power  in  the 
second-harmonic  mode  is 


dP2 

dz 


^wWImQr  dxj 

=  Wwd  Im(EEE  A*AfBm  exp[i(9"  +  0"  -  0*“)z  ] 

X  f  [<?„(*)  ]2<?j|"(x)  dx  \ ,  .  (11) 

■/  .u 


where  W  is  the  waveguide  width  in  the  y  direction. 
Equation  (11)  shows  immediately  that  a  phase-matched 
interaction  is  due  to  triplets  ( n,l,m )  or  space  harmonics 
for  which  the  exponent  in  Eq.  (11)  vanishes,  or,  using 
Eq.  (3),  a  phase-matched  interaction  occurs  when 

20"  -  &lu  +  (n  + 1  -  in){ 2tt/A)  =  0.  (12) 

Phase  matching  in  first  order  occurs  if  we  choose  the 
period  A  so  that  20"  -  0="  ±2ir/A  =  0.  The  synchronous 
contributions,  in  this  case,  arise  from  the  triplets 
(0,0,1),  (0 ,  —  1,0),  and  (-1,0,0).  A  simple  manipula¬ 
tion  of  Eq.  (11)  in  which  all  non-phase-matched  contri¬ 
butions  are  ignored,  in  which  0  is  assumed  to  be  equal 
to  the  bulk  propagation  vector,  and  where  the  energy 
is  assumed  to  be  confined  mostly  within  the  height  t  of 
the  guide  gives 


in  the  nondepleted  pump  approximation.  This  result  is 
of  a  form  identical  to  the  bulk  interaction6  except  that 
here  the  effective  nonlinear  coefficient  is 

dM»dJt(MX).  d«) 

The  approximate  equality  is  used  since  the  exact  nu¬ 
merical  coefficient,  which  is  of  the  order  of  magnitude 
of  unity,  involves  the  transverse  overlap  integral  and 
the  relative  magnitude  and  sign  of  the  three  first-order 
contributions.  It  will  be  included  in  a  more  complete 
forthcoming  paper. 

The  conversion  efficiency  from  w  to  2w  is  seen  to  be 
proportional  to  the  mode  power  density  Pu/Wt.  Since  W 
and  t  can  be  made  comparable  to  X,  this  power  density 
can  become  very  large  even  for  small  power  input.  The 
penalty  for  using  surface  corrugation  is  a  reduction  of 
the  effective  nonlinear  coefficient  by  the  factor  ~<71(.1/J’). 
Using  Eq.  (8)  in  the  case  of  a  GaAs  epitaxial  film  with 
/ rz  6  JU,  Wj=3.5,  and  «,-n3  =  0.2,  for  a  fundamental 
wave  at  \,  =  10.6  p  we  find  that  first-order  phase 
matching  requires  a  corrugation  with  a  period  of  A  =  107 
ju  and  that,  for  a  =  0.5  p,  drtt~^dCl Aa.  It  is  interesting 
to  note  that,  even  allowing  for  the  25 x  reduction  of 
the  effective  nonlinearity,  using  dQ^^l^xlO"21  MKS, 
the  effective  coefficient  is  comparable  to  that  of  LiNb03, 
one  of  the  best  phase-matchable  materials  presently 
used. 

Note  added  in  proof :  Nonlinear  interactions  involving  a 
spatial  periodic  modulation  of  the  nonlinear  coefficient 
d  rather  than  the  modulation  of  the  height  are  also  pos¬ 
sible  and  will  be  discussed  separately. 

Experimental  techniques  for  fabricating  surface  cor¬ 
rugations  in  GaAs  with  periods  as  small  as  0.28  p  have 
been  developed  in  cooperation  with  Dr.  H.  Garvin  of 
the  Hughes  Research  Laboratories.  The  availability  of 
such  techniques  plus  the  fast  evolving  technology  of 
GaAlAs  epitaxy  should  make  possible  the  development 
of  tunable  optical  parametric  oscillators,  unconverters, 
and  second-harmonic  generators  using  th’n  films. 
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The  Dower  threshold  tor  optieally  indueed  bulk  damage  in  sapphire  Is  a  sensitive  function  of 
the  ellipticity  of  the  incident  beam  shape.  Experimental  results  are  consistent  with  a  simple 

self-focusing  theory. 


We  wish  to  report  observations  of  the  strong  influence 
of  departures  from  circularity  of  an  optical  beam  on  the 
power  threshold  for  optically  induced  bulk  damage  in 
sapphire.  This  threshold,  which  is  known  to  arise  from 
th°  formation  of  a  catastrophic  self-focus  in  the  medi¬ 
um,1  is  found  to  increase  appreciably  when  the  beam 
cross  section  is  distorted  from  a  circle  to  an  ellipse, 
and  when  vertical  and  horizontal  confocal  parameters 
differ  (leading  to  ellipsoidal  phase  fronts).  These  ob¬ 
servations  are  consistent  with  the  theoretical  predic¬ 
tions  of  Vorob’yev2  and  Shvartsburg,3  which  are  extend¬ 
ed  slightly  here  to  apply  to  our  experiment. 

Our  experiments  were  performed  using  the  output  of  a 
single-mode  Q-switched  ruby  laser  and  amplifier  fo¬ 
cused  inside  sapphire  samples  with  different  lenses. 

The  laser  and  associated  monitoring  apparatus  are  de¬ 
scribed  in  detail  elsewhere.4  The  far-field  beam  pro¬ 
file  was  measured  to  be  Gaussian  down  to  8f7  of  the 
peak,  using  a  modified  multiple-lens  camera  technique. 
Typical  pulse  lengths  are  20  nsec.  The  sapphire  sam¬ 
ples  are  typically  3-in.  -long  by  i-in.  -square  bars. 

Damage  threshold  powers  for  circular  beams  of  differ¬ 
ent  sizes  are  shown  in  Table  I  and  compared  with  re¬ 
sults  for  elliptical  beams.  The  first  three  entries  in 
Table  I  show  the  increased  threshold  power  with  beam 
size  indicating  the  influence  of  electrostriction.1  The 
dimensions  listed  in  Table  I  are  those  at  the  beam 
waist,  near  which  the  self-focus  first  forms  at  thresh¬ 
old  for  circular  beams.  The  remaining  entries  in  Table 
I  show  the  effect  of  noncircular  beam  shapes.  The  fourth 
entry  gives  the  threshold  power  for  an  elliptical  beam 


whose  short  dimension  is  the  same  as  the  first  entry, 
whose  long  dimension  is  the  same  as  the  third  entry, 
and  whose  area  is  the  same  as  the  second  entry.  We  see 
that  the  threshold  power  is  appreciably  higher  for  the 
elliptical  beam  compared  with  any  of  the  circular 
beams.  The  last  two  entries  in  Table  I  show  data  for 
more  elongated  elliptical  beams  where  we  were  unable 
to  reach  threshold  at  the  maximum  power  available  from 
our  system. 

Each  power  threshold  in  Table  I  was  determined  from  a 
series  of  eight  laser  shots  of  differing  peak  powers. 
When  damage  occurred,  the  damaged  region  was  found 
to  possess  a  circular  rather  than  an  elliptical  cross 
section.  In  fact,  the  appearance  of  the  damage  tracks 
formed  above  threshold  was  indistinguishable  from 
those  caused  by  circular  beams.  The  tracks  themselves 
were  confined  to  the  region  between  the  two  line  foci  of 
our  astigmatic  optical  system,  but  there  is  no  evidence 
that  they  would  not  extend  beyond  the  upstream  focus  at 
higher  powers. 

To  correlate  these  observations  with  theory ,  we  have 
found  it  necessary  tc  resort  to  the  "paraxial-ray— con¬ 
stant-shape”  analysis  of  self-focusing  outlined  in  Ref.  5. 
More  accurate  numerical  solutions  of  the  nonlinear 
wave  equation,  such  as  those  employed  in  Ref.  1,  are 
extremely  difficult  to  obtain  in  this  case,  because  the 
noncircular  beam  shape  requires  an  additional  degree 
of  freedom  in  the  computer  code.  Nevertheless,  the  ap¬ 
proximate  analysis  has  had  some  qualitative  success, 
and  leads  to  simple  expressions . 
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TABLE  1.  Comparison  of  damage  threshold  powers  for  circular  and  elliptical  beams  of  different  sizes. 


Type  of  beam 

Radius1  (pm) 

Area  (pnr) 

Damage  threshold  power 
(MW) 

Astigmatic  focal  length11 
In  medium  (cm) 

Circular  beams 

14 

6.3  xio2 

0.  51  ±  0.04 

•  •  • 

(dimensions  at  low- 

'  37 

44  xio2 

0.  74  ±  0.  07 

•  •  • 

Intensity  waist) 

100 

308X102 

1.231  0. 10 

•  •  • 

Elliptical  beams  (dimen- 

14  X100 

44  xio2 

6. 0±  0.  50 

2.99 

slous  at  upstream  uw- 

14  X1200 

535  xio2 

>  17e 

OO 

Intensity  waist) 

7  X1200 

270  xio2 

>  17e 

OO 

‘The  beam  radius  is  defined  as  the  1/e  radius  for  the  intensity.  London,  1965),  p.  171]. 

'’This  is  the  distance  between  the  two  line  foci  uee  M.  Born  cUnable  to  reach  threshold  at  this  power, 

and  E.  Wolf,  Principles  of  Optics,  3rd  Ed.  (Pe  ‘gamon. 


Assuming  an  optical  field  of  the  form 

E  =  |£0  exp(i<p)  exp[i(ftz  -  wf)]  +  c .  c . , 

E0  =  Ejz)  exp[-  \(x*/a*  +  //ft2)  ] , 

we  find  by  the  method  of  Ref.  5  the  following  equations 
for  the  horizontal  and  vertical  beam  parameters  a  and 
ft : 


d2a 

1 

V 

(1) 

dz 2  a3 

~  ctb  ’ 

tf2ft 

1 

V 

(2) 

dz 2  ft3 

aft2 

Here  i]  =  P/Pl,  where  P  is  the  total  power  and  Pl  =n3c/ 
4eafea  in  Gaussian  units.  The  first  terms  on  the  right- 
hand  sides  of  Eqs.  (1)  and  (2)  represent  the  effect  of 
diffraction,  while  the  second  terms  arise  from  the  non¬ 
linear  contribution  ic^E3  to  the  dielectric  constant,  n  is 
the  linear  refractive  index.  Vorob’yev,2  who  derived 
Eqs.  (1)  and  (2)  in  a  different  way,  has  shown  that  they 
imply 

^(o2+fc2)=  0,  (3) 

which  can  be  integrated  to  give 

*V  +  ft2)  =  [feVo  +  ft2)  +  (2/oA)  (T]c  -  q)l z2 

+  2  /(a0a0  +  606a)z+A>2(a2+ft2).  (4) 

Here  a0  and  b0  are  the  axes  of  the  ellipse  formed  by  the 
e'1  points  of  the  transverse  intensity  profile  at  z=0, 
and  i]c  =  a0/2b0  +  ba/2a0.  The  initial  rates  of  change  a0 
and  b0  are  simply  related  to  the  principal  radii  of  curva¬ 
ture  Ra  and  Rb  of  the  phase  fronts  at  z  —  0: 

^«  =  ®o/Aj>  ■^»  =  ft(>/ftcr 

Using  Eq.  (4),  the  reader  may  easily  find  that  the  on- 
axis  intensity,  which  is  inversely  proportional  to  2 aft 
So2  +ft2,  becomes  infinite  at  the  point 

zfm=ka0b0(ri/Tic-l)'U2,  (5) 

for  plane  incident  phase  fronts  and  tj>  rjc.  It  is  therefore 
clear  that  the  critical  power  for  self-focusing  is 

Pc=nA  (*a =*»=«),  (6) 

which  equals  P1  for  circular  beams .  For  finite  values  of 
Ra  and  Rb,  a  real  solution  for  z  of  the  condition  a2  +b2 


=  0  occurs  only  if  P/Pl  exceeds 


For  P  =  t]ccR i,  the  self-focus  occurs  at  z=z0,  where 

-((PQ+bl)/z0  =  al/Ra  +  bi/Rb.  (8) 

For  higher  powers  the  self-focus  will  be  found  at 

z,~  1  +  (z0/ka0bQ)[(n  -  OAiJ1  n  ‘ 

For  the  fourth  entry  in  Table  I,  the  beam  parameters 
area0=0.179  mm,  fto  =  0.232  mm,  Ra  —  4.05  cm,  and 
R„  =  7.04  cm.  The  measured  power  threshold  for  the 
onset  of  bulk  damage  was  from  five  to  twelve  times 
higher  than  that  for  circular  beams  of  comparable  di¬ 
mensions,  and  ihe  damage  first  appeared  at  z  =  5. 8  cm 
from  the  entrance  face.  *  ;suming  that  the  damage 
should  first  appear  when  a  catastrophic  self-focus  forms 
at  P  =  VecPl,  we  expect  from  Eqs.  (7)  and  (8)  that  the 
threshold  power  should  be  increased  by  tjcc  =  12.7  and 
that  the  damage  at  threshold  should  appear  at  z0  =  5.51 
cm.  Thf  latter  figure  agrees  well  with  the  measured 
value,  and  the  change  in  critical  power  is  certainly  of 
the  right  order.  Unfortunately,  our  knowledge  of  the 
influence  of  the  transient  electrostrictive  contribution 
to  e2  is  too  poor  to  allow  us  to  know  which  of  the  first 
three  entries  in  Table  I  we  should  choose  as  Px.  The 
first  entry  (P^O.51  MW)  gives  excellent  agreement 
with  theory,  but  the  electrostrictive  response  for  an 
elliptical  beam  can  hardly  be  as  large  as  that  for  a 
circular  beam  whose  diameter  equals  the  small  radius 
of  the  ellipse.  Nevertheless,  the  agreement  with  theory 
is  sufficiently  good  to  allow  us  to  conclude  that  the 
damage  morphology  and  threshold  are  well  understood 
in  terms  of  the  self-focusing  theory. 

The  theory  outlined  above  gives  simple  formulas  for  the 
powers  and  distances  at  which  the  quantity  a2  +  ft2  van¬ 
ishes,  but  the  intensity  can  become  infinite  only  if  the 
area  nab  of  the  beam  van.  shes.  A  closer  examination  of 
the  solutions  of  Eqs.  (1)  and  (2)  shows  that  a  and  ft  both 
vanish  at  the  same  point.  Prior  to  this  point,  the  ratio 
a/6  approaches  unity.  This  is  consistent  with  the  ob¬ 
servation  that  the  cross  sections  of  the  damaged  region 
have  a  circular  shape  along  the  en*ire  length  of  the 
damage  track.  We  remind  the  reader  that  at  some  point 
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between  the  astigmatic  line  foci  the  low-intensity -beam 
cross  section  must  be  circular,  but  elsewhere  it  is 
elliptical 

These  results  imply  that  difficulties  which  arise  from 
bulk  damage  due  to  self-focusing  in  high-power  optical 
systems  can  be  avoided  or  ameliorated  by  the  proper 
choice  of  beam  parameters. 

The  authors  are  grateful  to  V.  Evtuhov  for  helpful  dis¬ 
cussions  and  for  providing  English  translations  of  Refs. 
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